Onset of DNA Aggregation in Presence of Monovalent and Multivalent 
Counterions ^ 

Yoram BurakQ Gil Ariel, and David Andelman 
School of Physics and Astronomy, 

Raymond and Beverly Sackler Faculty of Exact Sciences 
Tel Aviv University, Tel Aviv 69978, Israel 

We address theoretically aggregation of DNA segments by multivalent polyamines such as spermine 
and spermidine. In experiments, the aggregation occurs above a certain threshold concentration 
of multivalent ions. We demonstrate that the dependence of this threshold on the concentration 
of DNA has a simple form. When the DNA concentration cdna is smaller than the monovalent 
salt concentration, the threshold multivalent ion concentration depends linearly on cjjujAi having 
the form acoNA + /3- The coefficients a and /3 are related to the density profile of multivalent 
counterions around isolated DNA chains, at the onset of their aggregation. This analysis agrees 
extremely well with recent detailed measurements on DNA aggregation in the presence of spermine. 
From the fit to the experimental data, the number of condensed multivalent counterions per DNA 
chain can be deduced. A few other conclusions can then be reached: i) the number of condensed 
spermine ions at the onset of aggregation decreases with the addition of monovalent salt; n) 
the Poisson-Boltzmann theory over-estimates the number of condensed multivalent ions at high 
monovalent salt concentrations; in) our analysis of the data indicates that the DNA charge is not 
over-compensated by spermine at the onset of aggregation. 



INTRODUCTION 

Condensation and aggregation of DNA, induced by 
multivalent counterions, have been exte nsively studied in 
the pa st two decades (for a review, see iBloomfield et al.l 
l|200Cl|) and references therein). The term condensation 
usually refers to the collapse of a single, long DNA 
chain. Condensation plays an important role in stor- 
age and packing of D NA; for example, in viral capsids 
ijGelbart et all 1200(1 . Aggregation of DNA is a closely 
related phenomenon, where multiple chains attract each 
other and form a variety of condensed mesophases of com- 
plex structure l|Pelta et al.L llQQBalbT l . In both phenom- 
ena multivalent counterions play a crucial role, screen- 
ing the electrostatic repulsion between charged strands 
of DNA and mediating an effective attraction. 

A variety of tri- and tetra-valent ions can induce ag- 
gregation and condensation, among them the polyam ines 
sperm i dine (3"*") and spermi n e (4+ ) llChattorai et al. , 
1978[ iGosule and SchellmanL Il978t iTabor a nd Taboi . 
1984 ), as we ll as cobalt-hexamine llWidom and Baldwiii , 
198C , Il983^. In t ypical ex periments on ag grega- 
tion llPelta et all Il996bt iRasoaud etHl Il998t 
ISaminathan et all Il999(l multivalent ions are grad- 
ually added to a solution with fixed concentration 
of DNA segments and monovalent salt. Two such 
examples for spermine and s permidine are reproduced 
in Fig. 1 ijPelta et all Il996b|) . As the multivalent ion 
concentration is raised above a certain threshold, DNA 
segments begin to aggregate, and precipitate from the 
solution. Above the aggregation threshold, the DNA 
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FIG. 1 Percent of solubilized DNA, as function of polyamine 
concentration: square symbols - spermine, circles - spermi- 
dine. The solid and dashed lines are guides for the eye. DNA 
and NaCl concentrations are 3mM and 25 mM, respectively. 
Below the aggregation threshold, Caggr, and above the re- 
dissolution threshold, Credissoij_all the DN A is dissolved. The 
data is adapted from lPelta et al.l lll996bl) . 



concentration decreases gradually or abruptly, depend- 
ing on various parameters such as the monovalent salt 
concentration and total DNA concentration. Further 
addition of multivalent ions at higher concentrations 
reverses the aggregation. Above a second, re-dissolution 
threshold, all the DNA is re-dissolved in the solution 
(Fig. 1). The re-dissolution threshold (above which all 
the DNA re-dissolves) is almost independent on the 
DNA concentration. Its value can be attributed to 
screening of electrostatic interactions by multivalent ions 
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FIG. 2 Spermine concentration c^^aggr at the onset aggre- 
gation, as a function of DNA monomer concentration cdna- 
Data is shown for four monovalent sah concentrations: 2 mM 
(o), 13 mM (A), 23 mM (V), and 88 mM (□). The solid line 
corresponds to a fixed ratio: Cz . aggr/c pNA = 0.20. The data 
is adapted from lRaspaud etal] il998t) . 

l|Raspaud et al.Ul998|) . 

The aggregation threshold, where the onset of aggre- 
gation occurs, is the main experimental phenomenon ad- 
dressed in our theoretical paper. The multivalent ion con- 
centration at the onset depends strongly on the monova- 
lent salt and DNA concentrations. This dependence has 
been recently measured in detail for shor t (150 base pair) 
DNA segments in presence of spermine l|Raspaud et all 
Il998|) . and is reproduced in Fig. 2. The figure shows 
measurements of spermine concentrations at the onset of 
aggregation, for DNA concentrations ranging over four 
orders of magnitude and for four different monovalent 
salt concentrations: 2, 13, 23 and 88 mM. At very low 
DNA concentration, the spermine concentration depends 
strongly on the monovalent salt concentration. At higher 
DNA concentration it has only a weak dependence on the 
monovalent ion concentration but the spermine concen- 
tration is proportional to the DNA concentration, indi- 
cating that a certain number of spermine counterions are 
required, per DNA base, in order to indu ce aggregation. 
The s olid line in Fig. 2, adapted from iRaspaud et al.l 
|)1998D . corresponds to a ratio: c^.aggr/coNA = 0.20, 
where c^^aggr is the spermine concentration at the aggre- 
gation onset and cdna is the DNA concentration. This 
linear relation fits a large number of the experimental 
points in the intermed iate DNA conce ntration range. It 
has been suggested bv lRasoaud et al.l (^999, 1998) that 
the deviations from this line, at low and high DNA con- 
centrations, represent two distinct physical regimes that 
need to be analyzed separately from the intermediate 
regime, where the linear fit works well. 

In this work we focus on the onset of aggregation, and 
specifically on its dependence on the DNA concentration. 
We show that this dependence is simple for all the range 



of DNA concentration. Furthermore, for Cdna smaller 
than the monovalent salt concentration we show that this 
dependence is linear: c^^aggr = ckcdna+Z?- The coefficient 
(5 is the multivalent counterion concentration far away 
from the DNA chains, while a accounts for the excess 
of multivalent ions around each chain. These quantities 
can be extracted, e.g., from the four experimental curves 
of Fig. 2. Several further conclusions are then drawn 
on the onset of DNA aggregation and on the counterion 
distribution around each double-stranded DNA. 



THEORETICAL CONSIDERATIONS 

Consider an aqueous solution containing monovalent 
(1:1) salt, multivalent (z-.X) salt and DNA segments below 
their threshold for aggregation. Throughout this paper, 
the DNA solution is assumed to be dilute enough such 
that the DNA segments do not overlap. We also assume 
that these DNA segments can be regarded as rigid rods. 
The concentrations of added monovalent salt, multivalent 
salt and DNA monomers are denoted by Cs, Cz and cdna, 
respectively. These are the solute concentrations per unit 
volume as controlled and adjusted in experiments. We 
will assume that the monovalent and multivalent salts 
have the same type of co-ion, so that altogether there 
are three ion species in the solution: 

1. A multivalent counterion contributed from the z:l 
multivalent salt, of concentration c^. 

2. A monovalent counterion contributed by monova- 
lent salt of concentration c^, and by counterions 
dissociated from the DNA, of concentration cqna^ 
in total, Cs + cdna ■ 

3. Co-ions coming from both z:\ and 1:1 salts, of con- 
centration Cs + ZCz. 

Each DNA segment attracts a layer of oppositely 
charged counterions referred to as the condensed counte- 
rions. As long as the typical distance between segments 
is large compared to the electrostatic screening length 
K"^, the electrostatic potential decays exponentially to 
zero far away from the DNA segments. In turn, the con- 
centrations of the three ion species decay to well defined 
hulk values denoted by for the monovalent ions and 
for the z-valent ones. These concentrations should 
be distinguished from the concentrations Cs and Cz in- 
troduced above, which are the average concentrations of 
added salts regulated experimentally. 

The Debye screening length, characterizing the 

exponential decay of the electrostatic potential, is de- 
termined by the bulk concentrations of all three ionic 
species: 

^2 = 4^Zb [c^ -f z'cf + (c?° + zcD] . (1) 

where the third term is the co-ion concentration. It is 
equal to + zc^ due to charge neutrality far from 
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FIG. 3 Schematic representation of the muhivalent density 
profile, nz(f) between two neighboring DNA segments, each 
modeled as a cylinder of radius d. Here r is the distance from 
the axis of the left DNA strand. The radius r — R corresponds 
to the inter-strand mid-distance and is the unit cell radius. 
The density decays to its bulk value cf on distances larger 
than , where is the Debye length defined in Eq. Q 
The excess density of multivalent ions pz is indicated by the 
shaded areas. 



the DNA where the potential decays to zero. The 
above equation makes use of the Bjerrum length, Is — 

/ {ekBT), equal to about 7 A in aqueous solution at 
room temperature, fcsT is the thermal energy , e is the 
electron charge and e = 80 is the dielectric constant 
of water. The Debye length as well as are shown 
schematically in Fig. 3. Other quantities that will be 
defined below are also indicated in this figure. 

In dilute solutions different DNA segments do not over- 
lap. Following previous works, we introduce a cell model 
also shown schematically in Fig. 3. Note that the model 
serves to illustrate the subsequent derivations but is not 
essential for the validity of our main results. In the cell 
model, each segment, of a cylindrical cross-section, is 
at the center of a cylindrical cell of radius R and area 
A = ttE? such that 



CDNA = l/{aA). 



(2) 



Namely, each DNA monomer occupies a specific volume 
a A, where a ~ 1.7A is the average charge separation on 
the chain taken hereafter as the monomer length. 

We will assume below that the DNA solution is dilute 
enough so that R is large compared to the Debye length 
. This assumption is essential for our derivation and 
can be verified for all the experimental data considered 
in this paper. Density profiles of the three ion species are 
then practically identical to those near an isolated DNA 
segment with the same bulk concentrations c^j", c^. In 
other words, the profiles are determined uniquely by 
and c^, with practically no dependence (or, more pre- 
cisely, an exponentially small dependence) on the DNA 
monomer concentration. A demonstration of this claim 
is presented in Fig. 4, using the Poisson-Boltzmann the- 
ory in a cell model. For two very different values of R 
corresponding to different cdna, the counterion profiles 
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FIG. 4 Density profile nz(r) of 4-valent ions as function of 
r, the distance from the DNA axis, on a semi-log plot, calcu- 
lated using the Poisson-Boltzmann equation in a cell model, 
where the DNA segment is modeled as a uniformly charged 
cylinder. The cell model is shown schematically in the in- 
set. Two cell sizes are shown, with outer radii Ri = 560 A 
(cdna = ImM) and R2 = 1.8 x 10* A (cdna = 10"^ mM), 
indicated by arrows. In both cases the radius of closest ap- 
proach of ions to the charged chain is at r = d, where d = 10 A 
as indicated by a dotted vertical line. The boundary condi- 
tion at the inner cylinder matches the linear charge density of 
DNA (le/1.7A). The bulk densities of monovalent and mul- 
tivalent ions, cj° and c^, are chosen to be the same in the 
two cells, leading to practically identical density profiles. The 
solid line represents the larger cell {R2), and diamonds are 
used for the smaller cell (-Ri). Density profiles of monovalent 
counterions and co-ions are not shown but are also practi- 
cally identical in the two cells. Average salt concentrations 
are Cs = 22 mM and Cz = 0.21 mM in the smaller cell, and 
Cs = 23 mM, Cz = 0.039 mM in the larger cell. Bulk con- 
centrations are cf = 23 mM and cf = 0.039 mM. Note that 
these bulk concentrations are practically identical to the salt 
concentrations in the larger cell. Note also that c^" > Cs in 
the smaller cell reflecting the contribution of the counterions 
released by the DNA. 



match perfectly when the values of and are the 
same. Note that the average concentrations of added 
salts, Cs and Cz, have different values in the two cells 
because of the contribution of condensed ions. 

The total number of z-valent counterions, per cell unit 
length, is given by: 



Ac, = Ac° 



(3) 



where is the excess number of z-valent ions per unit 
length near the DNA. Throughout the paper we use the 
symbol c to denote concentrations per unit volume and 
p for concentrations per DNA unit length. The excess 
Pz can be evaluated in the limit of infinite cell radius, 
corresponding to an isolated chain: 



Pi 



rdr [uz (r) — c° 



(4) 
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where Uzir) is the z-valent local counterion concentration 
at distance r from the axis of symmetry, and Uzioo) — 
d^. Following the discussion in the previous paragraph, 
the excess pz is determined uniquely by c^j" and cf. Its 
exact functional dependence on these variables is gener- 
ally not known, although it can be evaluated approxi- 
mately, e.g., using the Poisson-Boltzmann equation or in 
computer simulations. 

For monovalent counterions we have, in a similar fash- 
ion: 

Ac,+AcuNA = Acr+PiicT,cT), (5) 

where pi , the excess of monovalent counterions per unit 
length, is defined as in Eq. 0] and Acdna = 1/a is the 
DNA charge density per unit length. The extra term in 
the left-hand-side of Eq. |31 originates from monovalent 
counterions contributed by the DNA monomers. Using 
Eq. [51we can rewrite Eqs. |31and|5las: 

Cz=cf + apz {cf , Cf ) CDNA (6) 

and 

Cs=c^ + [api (cj°, c^) - 1] CDNA. (7) 

These two equations relate the experimentally adjustable 
Cs , Cz and cdna to the bulk densities , that in turn, 
are important because they determine the ion density 
profiles. 

In the limit of infinite DNA dilution, cdna = 0, and 
therefore Cz = and Cg = c^. At any finite DNA 
concentration Cz and Cg are not equal to and cj°, 
respectively, because each segment captures some of the 
multivalent ions and releases a number of monovalent 
ones. Equations|Sland[7|express the correction to Cg, Cz at 
given c^, for both mono- and multi-valent counterion 
species. The dimensionless quantities api,apz are the 
excess of the mono- and multi-valent counterion species, 
respectively, per DNA monomer. 

We would like to emphasize the generality of Eqs. El 
and[7| They do not depend on the assumption of parallel 
DNA residing in the middle of oriented cylindrical unit 
cells, or on any mean-field approximation for the distri- 
bution of counterions. The only assumption required to 
derive Eqs. Eland[7|is that the average distance between 
DNA segments is large compared with the Debye length. 
Although Eqs. and [7| are correct for any Cg, Cz and 
Cdna below the onset of DNA aggregation, we will be 
interested below specifically in the aggregation onset. 



Onset of aggregation 

Our aim now is to find how the value of Cz at the onset 
of aggregation, Cz, aggr, depends on cdna- We will assume 
that this aggregation onset depends on and c^, but 
not on the average distance between DNA chains. We 
motivate this assumption by the fact that c^" and 



determine the density profile of multivalent counterions 
around the DNA chains, which, in turn, mediate the at- 
traction necessary for aggregation. Before discussing this 
assumption in more detail, let us first consider its con- 
sequences. We can imagine an experiment where is 
gradually increased while c^j" is kept fixed. Aggregation 
will start, in this experiment, above a certain threshold 
value of c^. Our assumption is that this threshold does 
not depend on cdna- In real experiments, however, Cz is 
adjusted rather than c^, and Cg is kept fixed rather than 
c^. In order to find the threshold value in terms of the 
experimentally available Cz we need to map c^, onto 
Cg,Cz. This mapping is described by Eqs. |S]-|7| and in- 
volves cdna- It is only through this mapping that cdna 
will affect the threshold of aggregation. 



The limit of cdna ^ Cg : 

The limit cdna ^ Cg offers a particularly simple de- 
pendence of Cz^aggr on CDNA and is considered first. Most 
models and experiments indicate that monovalent coun- 
terions cannot overcharge DNA segments. Hence the 
monovalent excess, api, in Eq. 13 is a number between 
zero and one, because the excess monovalent charge is 
smaller than that of DNA. From Eq.[71|cs — <C Cs as 
long as Cdna <^ c^. It is then possible to replace c^ by 
Cg, leading to a simplification of Eq. |^ 

Cz = c^ -I- apz (cs, cf') Cdna- (8) 

Note that cdna is indeed smaller than Cg in most of the 
experimental points in Fig. 2. However a similar simpli- 
fication cannot be applied for Cz because it is typically 
much smaller than Cg, and often smaller than cdna- 

According to our principal assumption, aggregation 
starts at a threshold value c^ = c*, which does not 
depend on cdna (while Cz. aggr, the average multivalent 
salt concentration does depend on cdna through Eq.|SJ|. 
Similarly, the density profile at the threshold does not de- 
pend on Cdna, because it is determined by c"^ — Cg and 
c*. The excess of z-valent counterions, as determined 
from this profile, is equal to: 

P*z= pAcs,c*z), (9) 

with no dependence on cdna- Using the threshold values 
c* and Pz in Eq.|Sl we find that the average concentration 
of z-valent ions at the onset of aggregation is: 

Cz, aggr(cDNA) = C* + aPzCDNA- (10) 

This is the threshold concentration that was measured 
experimentally in RasDaud ct al. (1998). Note that in 
Eq. ^1 c* as well as p* depend on the monovalent salt 
concentration, Cg, but the explicit dependence is omitted 
for clarity. 

The simple relationship expressed by Eq. 1101 is one of 
our main results. As a visualization of this result we refer 
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FIG. 5 Spermine concentration at the onset of aggregation c^^aggr as a function of cdna, fitted to the form derived in Eg. IIUI 
(different line types a re used for different sal t concentrations) . Value of Cs (in mM) is indicated next to each curve. Experimental 
data is adapted from iR.aspaud et alJ (jlQ^ST) and shown in the following symbols: Cs = 2mM (o), 13 mM (A), 23 mM (v), and 
88 mM (□). Experimental error bars (E. Raspaud, private communication) are indicated by vertical lines. The fitted lines and 
experimental points are shown using a linear scale in (a) up to cdna = 1-5 mM, and a log-log scale in (b) up to cdna = 100 mM, 
allowing all data points to be shown on the same plot. Only the data up to cdna = 10 mM was used for the linear fit. The 
crossover values of cdna, as defined by Eg. 1141 are indicated by arrows in (b). 



again to Fig. 3. The quantities pz, cf and the density 
profile nz{r) are indicated in this figure. At the onset 
of aggregation is equal to c* and does not depend 
on Cdna (or equivalently, on the spacing between DNA 
segments, R). As cdna is increased the distance between 
DNA strands decreases. The onset values of and 
do not change, but the contribution of pz to the average 
concentration increases, leading to an increase in Cz^aggr- 
The coefficients op* and c* of the linear dependence 
in Eq. ^| are the coefficients a and /9 defined in the in- 
troduction section. They can be easily found from the 
experimental data: c* is the value of c^^aggr in the limit 
of infinite DNA dilution, cdna — >■ 0, since in this limit 
Cz = = c*. The excess at the onset, p*, can be 
found from the slope of Cz,aggr as function of cdna- Be- 
fore presenting a detailed comparison with experiments, 
we generalize the treatment for small cdna to arbitrary 
values. 



The case of cdna > Cg : 

When Cdna is of the same order as Cs or larger, correc- 
tions to c'^ must be taken into account, as expressed by 
Eq. [71 and the linear relation of Eq. ^|no longer holds. 
The ion density profiles as well as Cs and Cz are now de- 
termined by the two variables and c^. The relation 
between c'^ and and the experimentally controlled 
Cs, Cz, Cdna is given by Eas. l6l7l In terms of c^j", c^ the 
criterion for aggregation remains the same as in the pre- 
vious case: 

cT = <{cT)- (11) 



The three equations IHl ITIand ll II with the three unknowns 
, and Cz lead to a unique solution for Cz, aggr- Note 
that is larger than Cs because of counterions coming 
from the DNA as can be seen in Eq. [7| where api — 1 is 
negative. In Eg. IIUI is replaced by c^j", which is larger 
than Cs for large cdna- Hence, increasing cdna has an 
effect similar to addition of monovalent salt. As noted 
above, this effect is significant only for cdna > Cg. 



COMPARISON WITH EXPERIMENT 

iRaspau d et al.l (^9^ measured the spermine {z = 4) 
concentration Cz at the onset of aggregation as a function 
of Cdna for four values of Cg and with cdna ranging over 
four orders of magnitude — from 10~^ to 10^ mM. We 
fitted the data (E. Raspaud and J.-L. Sikorav, private 
communication) for each Cs to a straight line according 
to Eq. ^1 The least square fit presented in Fig. 5 takes 
into account the experimental error bars and data points 
up to Cdna — 10 mM. In Fig. 5 a the fit is shown using 
a linear scale which covers the range of cdna only up to 
Cdna = 1.5 mM for clarity purposes. Due to the large 
range of cdna it is impossible to show all the data on the 
linear scale of Fig. 5 a. Instead, the same data and linear 
lines are shown in Fig. 5 6 on a log-log scale over the full 
experimental range of cdna- 

The linear fit is very good for all four values of monova- 
lent salt concentration Cg . Note that for Cg — 88 mM the 
fit is very good up to the largest value of cdna — 48 mM 
reported in the experiment, although our fit takes into 
account only data poi nts up to cdna = 1 mM. It was 
previously suggested (jRaspaud et al.l Il998j) that a sepa- 
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rate regime exists for cdna ^ 10 mM, characterized by a 
power law relation between and cdna with an expo- 
nent smaller than unity. Our analysis suggests a different 
conclusion. The fit clearly demonstrates that the relation 
is linear all the way up to cdna — 48 mM, as predicted 
by Eq. EH Note also that even at cdna = 48 mM we 
have Cdna < Cs so the assumptions leading to Ea. llUI are 
still valid. 

The only points in Fig. 5 h that deviate significantly 
from the fit are the three points where Cg = 13 mM (tri- 
angles) and Cdna > 20 mM (two of these points coincide 
with points having Cg — 88 mM, shown using square sym- 
bols.) This deviation is easily explained by the fact that 
Cdna S> Cs so that corrections to must be taken into 
account. For example, at cdna — 90 mM the nominal 
monovalent counterion concentration is 103 mM, taking 
into account counterions contributed by the DNA. In or- 
der to find cj° we need to subtract the condensed coun- 
terions, as determined by p\. We can estimate p\ at 
this point by solving the Poisson-Boltzmann equation in 
a unit cell with the appropriate radius. The chemical 
potentials of the three ion species are tuned such that 
their concentrations match the known values of and 
Cs- This leads to an estimate: ~ 68 mM. Hence, 
Cz at the onset of aggregation should lie a little below 
the continuation of the Cs = 88 mM line which is, in- 
deed, where it is found. The trend for = 13 mM can 
probably be seen already at the point cdna — 15 mM, 
although the deviation at this point is still within the 
range of experimental error. The few other experimental 
points with cdna ~ Cs deviate slightly from the straight 
line as well (still within experimental error bars). In all 
these cases the deviation is in the direction corresponding 
to a higher value of c^., as expected. 

A linear relation of the form c^^aggr = ckcdna + P, was 
previously suggested on empirical basis for aggregation 
induced by spermidine (3"*"), on a smaller range of DNA 
concen trations (Osland and Kleppe . 1972.; Pelta et al.L 
ll996tJ) . Although this result looks similar to our predic- 
tion on the onset of aggregation, it is not directly related 
to our analysis because c^^aggr was taken in those works 
to be the transition midpoint. This is the point where 
half of the maximal precipitation of DNA is reached. Our 
analysis does not apply at the transition midpoint since it 
requires all the DNA segments to be well separated from 
each other. Indeed, the coeffi cient ct, related to the tran- 
sitio n midpoint, was fou nd in lOsland and Klepod l)l977|) 
and 'Pclta et alJ l)l996bf) to be of order 10^, much larger 
than unity. Such a value of a cannot be interpreted as 
the excess of spermidine ions per monomer near isolated 
chains. 

The parameters of the linear fit in Fig. 5 are summa- 
rized in Table 1 for the four experimentally used values 

of Cs. 



TABLE 1 Fit parameters used in Fig. 5. 



Cs [mM] 


cl[mM] 




2 


± 0.0003 


0.194 ±0.020 


13 


0.011 ±0.002 


0.191 ±0.013 


23 


0.031 ± 0.005 


0.173 ±0.025 


88 


0.52 ±0.05 


0.135 ±0.026 



Crossover in the log-log plot 

For presentation purposes we plot in Fi g 5 &, c^ ay^r vs. 
cpNA on a log-log scale, as appeared in iRaspaud et alJ 
l)l998|) . The linear relation that was found between these 
two quantities is not clearly manifested on the log-log 
plot, because a linear dependence of the form y = ax + (3 
is not easily recognized in such a plot. Furthermore, such 
a linear relation appears on a log-log plot to be artificially 
characterized by two distinct behaviors, at low and high 
values of the inde pendent variable. The se two behaviors 
were mentioned in lRaspaud et al.1 l|l998|) and can be seen 
in Fig. 5 h. However, they do not represent in our opin- 
ion two real physical regimes and can be understood by 
taking the logarithm of Eq. EH For small cdna (large 
R): 

log -log c: (12) 

i.e, Cz does not depend on Cdna as is seen in Fig. 5 & in 
the small cdna limit. In the opposite limit of large cdna 
(small R): 

log Cz ~ log CDNA + log apl (13) 

Here, the linear dependence of Cz on cdna yields a line 
of slope 1 in the same figure. 

The crossover between these apparent behaviors occurs 
when the number of bulk and excess ions are the same: 



When Cdna is much smaller than this crossover value, the 
number of excess multivalent ions near DNA segments is 
negligible compared to their total number. In the other 
extreme of cdna much larger than the crossover value, 
the number of free multivalent ions is negligible compared 
to the excess ions, and nearly all multivalent ions are 
bound to the DNA. 

For the experimental data in Fig. 5 the crossover value 
is equal to 0.06, 0.18 and 3.9 mM for = 13, 23 and 
88 mM, respectively, and smaller than 1.5 x 10"'^ mM for 
Cs = 2 mM. The first three crossover points are indicated 
by arrows in Fig. 5 b. 

DNA AGGREGATION AND COUNTERION 
CONDENSATION 

We separate the discussion following our results in 
three parts. The first addresses the conditions required 
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for DNA aggregation. The coefRcients of the hnear rela- 
tion in Eq. ^1 c* and p*, have a definite physical mean- 
ing. Their values, as extracted from the experimental 
data provide insight on these conditions. The second 
part deals with condensation of countcrions on DNA (to 
be distinguished from condensation of DNA chains). The 
general relation pz = Pz{c'^,c'^) that was introduced in 
Eqs. 13 - 21 is a property of counterion condensation on 
isolated chains. By extracting the values of pz, and 
at the onset of DNA aggregation, we can learn about 
exact density profiles of spermine around DNA, and com- 
pare our findings with approximations such as Poisson- 
Boltzmann theory. Finally, we comment on our main 
assumption, which was used in the theoretical consider- 
ations section. 
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Conditions at the onset of aggregation 

Most of the proposed theoretical models for inter- 
chain attraction an d aggr eg ation (see, for exam- 
ple, lArenzon etll] lll999|: iBoriikhov et all ll2n0lL 
2n02l): iHa and Liul hmf\: 'Nguve n et all |2n0(| : 

Olvera de la Cruz et al.t (,1995): Raspaud et al.l l)l998|) : 
Wittmer et alJ l)l995f) ') regard the charged chain as 
surrounded by a layer of condensed ions which is usually 
modeled as a one-dimensional gas. This layer mediates 
an inter-chain attraction, and the models predict the 
number of condensed ions required to initiate aggrega- 
tion of the chains. In the current work we do not address 
this theoretical problem, but rather concentrate on what 
can be inferred from the experimental results using the 
analysis presented in the previous section. This analysis 
provides insight on the conditions prevailing at the onset 
of aggregation. In particular, the excess p* characterizes 
the number of condensed multivalent counterions that 
are present near each chain at the onset. Although in 
general the notion of condensed counterions is somewhat 
ill-defined, as it depends on which ions are regarded as 
bound to the DNA, we show in the Appendix that in 
our case it does have a reasonably well defined meaning. 
Furthermore, the number of condensed multivalent ions 
per monomer is practically the same as ap*. 

The excess of multivalent counterions per monomer, 
ap*, is shown in Fig. 6 as function of Cg. All values are 
taken from Table 1, as extracted from the experimental 
data. The dashed line is a linear fit. Two different axis 
scales are used on the left and right of the plot. The left 
axis shows the value of ap*. The right one shows the 
part of DNA charge that is compensated by condensed 
multivalent ions, zapl, where z = 4 for spermine. From 
the plot we deduce the following two conclusions: 

1. The number of condensed multivalent ions (per 
DNA monomer) ap* at the onset decreases as the 
monovalent salt concentration increases, with vari- 
ation between 0.19 and 0.14. A possible reason for 
this trend may be that the bare electrostatic repul- 
sion between chains is decreased due to increased 



FIG. 6 Excess of multivalent counterions per monomer at 
the onset of aggregation, apl, as function of Ca. All values 
are tak en from Table 1 , as ex tracted from the experimental 
data of iRaspaud et all lll998t) . Error bars are indicated by 
vertical bars and the dashed line is a linear fit to be used as 
a guide to the eye. On the right axis zap* is shown, where 
z — 4 for spermine. This value is equal to the fraction of DNA 
charge compensated by the condensed multivalent ions. Note 
that according to the Manning condensation theory the same 
quantity is equal to 0.94, for tetravalent ions and no added 
salt. 



screening. Hence a smaller number of multivalent 
ions is required in order to overcome this repulsion. 
The change in p* may also be related to the com- 
petition between monovalent and multivalent ions 
in the aggregated DNA state. 



2. The data indicates that there is no over-charging 
of the DNA by spe rmine at the onset (see also 
iNguven et alJ pOOOj) ') since zap* < 1. At higher 
concentration of spermine, beyond the threshold, 
we do not rule out the possi bility of D NA over- 
charging, as was suggested bv iNguven et al... (2000.') . 



Although p* decreases with increase of Cs, it is of 
the same order of magnitude for all the Cs values in 
Table 1. In contrast, c* varies in Table 1 over more 
than th ree orders of magnitude. As wa s previously sug- 
geste d l|01vera de la Cruz et al.l Il995l: iRasoaud et al.l 
Il998|) . this large variation in c* is a result of compe- 
tition between monovalent and multivalent counterions. 
We discuss the relation between p* and c* to some ex- 
tent in the following subsection. A more detailed analy- 
sis of this relation, emphasizing the role of competition 
between the two counterion speci es, will be presente d 
in a separate publication (see also. iBelloni et al.) lll984t) : 
IWilson and Bloomfieldl ^ll979^ : IWilson et all (119801) ). 
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TABLE 2 Excess of 4-valent ions near DNA 
compared with PB theory. 



cf=[mM] 




ap^ (exp) 


apz (PB) 


2 


± 0.0003 


0.194 ±0.020 


0.186 ±0.005 


13 


0.011 ±0.002 


0.191 ±0.013 


0.178 ±0.002 


23 


0.031 ± 0.005 


0.173 ±0.025 


0.172 ± 0.002 


88 


0.52 ±0.05 


0.135 ±0.026 


0.164 ±0.002 



Counterion condensation 

We now turn to analyze the condensation of monova- 
lent and multivalent ions around DNA. Each line in Ta- 
ble 1 provides a measurement of the excess pz at certain 
values of cf and cf. The general relation pz{c^,cf) is 
a property of counterion density profiles around isolated 
DNA segments. Hence, the data in Table 1 can be used 
to test any particular theory used to calculate such ion 
distributions. 

The most simple model to cons ider is the P oisson - 
Boltzmann (PB) th eor y (se e fAndelman (1994') ; 
iGueron and Weisbuch I l)l98(]|) : .Le Bret and Zimm 
l|l984j) : lOosawal l)l97l[l ). In Table 2 we compare the 
excess predicted by PB theory with the experimental 
result, by solving the PB equation such that cf and 
match the experimental values of and c* from Table 1. 
The excess is then calculated from the PB density 
profile, and compared with the experimental value of 
apz (equal to ap* of Table 1). The DNA is modeled as a 
uniformly charged cylinder of radius d = 10 A. 

Inspection of the results in Table 2 shows that there is a 
reasonable agreement with experiment (within the error 
bars) for the three smaller values of Cg = 2, 13, 23 mM. 
However, for Cg = 88 mM there is a 30% deviation. The 
two data points with cdna > 10 mM that were not taken 
into account in the linear fit of Fig. 5 suggest that pz is 
closer to the lower bound of the experimental error range, 
whereas the PB value is larger than the upper bound. 

Overall, the agreement with PB theory (Table 2) is 
surprisingly good considering that PB theory does not 
work so well for bulky multivalent ions. Deviations from 
PB theory have several sources. One of these sources is 
specific molecular details such as the geometrical shape of 
ions, DNA structure and short-range interactions. An- 
other source for deviations are ion-ion correlations be- 
tween spermine molecules, computed in theories which 
go beyond the mean-field approximation. However, these 
correlations tend t o increase the number of bound multi - 
valent counterions ijLvubartsev and Nordensk iold'.'l997*). 
while for Cg = 88 mM, the number of bound multivalent 
counterions is decreased. We conclude that correlation 
effects by themselves are not the main source of the de- 
viations seen in Table 2. In addition the data analy- 
sis does not indicate over-charging of the DNA. Such an 
effect may be expec ted if correlation effects are strong 
ijNguven et al.ll2000() . 

In Fig. 7 we compare the DNA aggregation data with 
PB predictions at finite DNA concentrations. For each 



DNA concentration the PB equation is solved in a cylin- 
drical cell of the appropriate radius. The multivalent 
counterion concentration Cz is gradually increased until 
the onset is reached, and its onset value, Cz^aggr is plotted 
as function of cdna- Two different criteria are used to 
determine the onset Cz^aggi - In Fig. 7 a it is chosen as the 
point where is equal to the experimental value c* of 
Table 1; whereas in Fig. 7 b the onset is chosen the point 
where pz — p*. In order to span all the data range we 
use for convenience a log- log plot, as in Fig. 5 b. 

On a linear scale all the lines in Fig. 7 a and b are 
straight lines. This fact serves as additional confirmation 
of our general analysis in the theoretical considerations 
section. In accordance with our analysis, both c* and pz 
are constant along each line, and the slope of each line is 
equal to apz ■ Note that the relation between c* and pz is 
determined in Fig. 7 within the PB approximation, while 
in Fig. 5 both of these coefficients are related to the actual 
counterion density profiles in the experimental system. 
The use of the PB equation is the source of deviations 
from experimental data in Fig. 7. 

On first inspection the match with experiment in 
Fig. 7 a is very good, whereas the match in Fig. 7 6 is 
not as good. On closer inspection it is seen that the fit in 
Fig. 7 6 is not good for small values of cdna, while it is 
actually better than in Fig. 7 a for large cdna- With the 
PB equation it is not possible to obtain a perfect fit for 
both small and large cdna because the values of Cz and 
Pz are not independent. Fixing = c* (as in Fig. 7 a) 
sets a value of pz that is different from p*; and the op- 
posite happens in Fig. 7 b. The fit in Fig. 7 a is quite 
good even for large cdna because the values of p* are of 
similar order of magnitude for all four lines. 

Deviations as in Fig. 7 are inevitable if any approx- 
imations are used to model the distribution of counte- 
rions around DNA. Note however that within such ap- 
proximate models our general theoretical considerations 
should apply, as long as the total number of ions in the 
system is counted proper ly. Such a model that goes be - 
yond PB was proposed in lNguven and Shklovski il ll2fl0l|^ . 
Indeed, within this model a linear relationship similar to 
Eq. ^1 was found. 

The experimental results analyzed in this section may 
be influenced, to a certain degree, by the fact that there 
was more than one type of monovalent counterion in the 
system. For the three higher salt concentrations, except 
for Cs — 2 mM, the solut ion contained 10 mM of TrisH"*" 
ions in addition to Na+ ijRaspaud et al.i ri998). For the 
largest salt concentration, 88 mM, where significant devi- 
ations from PB theory are found, this effect is probably 
negligible. Another detail regarding the TE buffer is that 
the Tris ions may be only partly ionized. If only 80% 
of Tris is ionized, as suggested in Tang ct al. (1997j), the 
concentrations Cg = 13 mM, 23 mM and 88 mM should be 
reduced by 2 mM. Although this will have only a small 
effect on our results, it will improve both the comparison 
with PB and the fit with the dashed line in Fig. 6, for 
the point Cg = 13 mM. For the two other concentrations 
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FIG. 7 Spermine concentration (in mM) as a function of DNA monomer concentration (mM) at the onset of aggregation, 
calculated using the PB equation. Two different criteria are used in parts (a) and (b) to determine the onset: in (a) c^, as 
calculated using the PB equation, is equal to the experimental value of from Table 1. In (b) pz of PB theory is equal to 
p* from Table 1. The radius of DNA is taken as d = 10 A. Log- log plot is used in order to show the five decades of DNA 

concentrations. For each Cs the plot covers experimental data up to cdna = Cs. For larger cdna, corrections due to changes in 
cj° should be taken into account, as was discussed in the preceding section. All notations are the same as in Fig. 5. 



of 23 mM and 88 mM the effect will be negligible. 



Further comments on underlying model assumption 

Our underlying assumption, that the onset of aggrega- 
tion depends uniquely on and (but not on Cdna), 
is an approximation that can be justified on several differ- 
ent levels but deserves further and more thorough investi- 
gation. The most simple motivation for this assumption 
is that multivalent ions, in the vicinity of the chains, me- 
diate the attraction necessary for aggregation. In turn, 
the number of condensed multivalent ions near each chain 
is determined by and c^. 

Let us first suppose that aggregation starts when a net 
attraction appears between two chains. This assumption 
may be justified if chains are sufficiently long and their 
translational entropy can be neglected. In order to find 
the onset of two-chain attraction the free energy of a two- 
chain complex should be calculated as a function of the 
distance between the two chains. This free energy rep- 
resents the effective interaction between the two chains, 
mediated by the ionic solution. The counterion distri- 
bution near each chain will not be the same for close-by 
and for isolated chains. However in both cases the con- 
centrations must decay to their bulk values throughout 
the solution, cj" and d^. This requirement serves as a 
boundary condition, imposed at a large distance from 
the two chains. It will determine uniquely the counte- 
rion distribution between the chains, as well as the free 
energy associated with the two-chain complex. Hence 
cj° and determine the effective interaction between 
chains, and in particular whether an attraction occurs 
at a certain range of inter-chain separations; in terms of 



these variables the onset of two-chain attraction does not 
depend on cdna- 

Strictly speaking, the onset of aggregation and the on- 
set of two-chain attraction are not the same. The aggre- 
gate phase involves interactions between multiple chains, 
whereas chains in the dilute phase interact very weakly 
with each other. Aggregation starts when the free en- 
ergy per chain is equal in the dilute and aggregate phases. 
Note that the chemical potential of each ion species must 
be the same in the two phases, and that in the dilute 
phase these chemical potentials arc directly related to 
and . Hence cf and determine the free energy 
per chain in the two phases. The approximation of inde- 
pendence on Cdna neglects the translational entropy of 
DNA segments, which can be justified for long enough 
and rigid segments. It also neglects contributions from 
interactions between chains in the dilute phase, which 
are assumed to be small compared to the free energy of 
the single DNA-counterion complexes. 



SUMMARY 

We have shown that the onset of aggregation at finite 
(non-zero) DNA concentration, Cz^aggr, is determined by 
the onset in the limit of infinite DNA dilution. For DNA 
monomer concentration smaller than that of monovalent 
salt, Cdna ^ Cs, the multivalent counterion concentration 
at the onset, Cz^aggr, depends linearly on Cdna- The coef- 
ficients of this linear dependence arc the bulk concentra- 
tion of multivalent counterions and their excess relative 
to the bulk near each DNA segment. Both of these co- 
efficients are of theoretical interest and can be extracted 
from the available experimental data. 



Our main assumption is that the onset of aggregation 
can be related to the ion density profiles around each 
chain. Hence, it is uniquely determined by cf and cf, 
the bulk concentrations of the two counterion species, 
respectively. Our results and fit to experiment strongly 
support this assumption. Nevertheless, we believe that 
more detailed theoretical and experimental investigations 
are needed in order to fully understand its range of va- 
lidity. For example, it will be of interest to test experi- 
mentally the equilibration of a DNA solution through a 
dialys is membrane, with a cell containing only counte- 
rions fBraunl in et al.l ll982t lPl um and BloomfjeldLlToSSt 
ISubirana and Vivesl Il98lj) This procedure allows a di- 
rect control of the ionic bulk concentrations. 

In order to predict precisely the onset of aggregation, 
the structure of the aggregated phase must be consid- 
ered. Nevertheless, it is instructive to focus only on single 
chains at the onset, as is often done. At the aggregation 
onset the electrostatic repulsion between isolated chains 
in solution must be overcome by a sufficiently strong at- 
traction mediated by multivalent counterions. This num- 
ber of counterions is expected to depend only weakly on 
physical parameters such as the monovalent salt concen- 
tration. Our analysis does not address directly the ques- 
tion of the onset origin, but merely supports the fact that 
the number of condensed multivalent ions at the onset, 
apl, is of the same order of magnitude, regardless of the 
Cs value. A more refined result of our analysis is that 
ap* is not constant but decreases with increase of Cg- On 
the other hand c*, the value of at the onset, depends 
strongly on Cg. This is mainly a result of the competi- 
tion between monovalent and multivalent ions, as will be 
addressed in a separate publication. 

Our analysis also sheds light on counterion condensa- 
tion on DNA, which is independent on the criterion for 
DNA aggregation. The experimental data indicates that 
for high Cs the number of spermine ions in the vicinity of 
DNA is smaller than the prediction of Poisson-Boltzmann 
theory. A similar trend was ob served in computer simu- 
lations (Lvubartsev and Nordenskioldl Il997t) of spermi- 
dine (S"*") and NaCl in contact with DNA. Spermidine 
binding was affected by addition of monovalent salt more 
strongly than the Poisson-Boltzmann prediction. For 
high salt concentrations spermidine binding was consid- 
erably smaller. In the computer simulations both molec- 
ular specific interactions, the geometrical shape of the 
constituents and inter-ion correlations were taken into 
account. All these effects, and in particular the geome- 
try of the spermidine molecule, which is similar to that 
of spermine, were found to play an important role. 

The above analysis demonstrates that specific in- 
teractions play an important role in determining the 
threshold of aggregation. In the dilute phase these in- 
teractions strongly influence the competition between 
monovalent and multivalent ions and the free energy 
of DNA-counterion complexes. Similarly, specific in- 
teractions p l ay a prominent role in the dense phase 
ijStrev et al.Ul99^ . Force measurements under osmotic 
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stress l)Rau et al 1 11984 iRau and Parseeianl Il992albt) 
provide a wealth of information on these interactions. 

In conclusion, the physical parameters extracted here 
from experiment on the onset of DNA aggregation pro- 
vide insight on the conditions required for aggregation, 
and on condensation of ions around DNA. These param- 
eters may turn out to be of great value in assessment 
of various theoretical models. Additional detailed exper- 
iments may further deepen our understanding of these 
complex phenomena. 
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APPENDIX 

In this appendix we discuss the relation between the 
excess and the number of condensed ions. The latter 
quantity is not as well defined as the former, but relates 
more naturally to the aggregation mechanism. The no- 
tion of condensed ions suggests that some ions are bound 
to the charged chain while others are free. In reality there 
is a density profile that extends all the way from r = d 
to r = R with no definite separation between condensed 
and free ions. In the following we define condensed ions 
rather loosely as the number of ions up to a certain char- 
acteristic dis t ance from the chain (^elloni et al., 198-^ 
I Wilson et al.L IT98(il) . We show that for multivalent ions 
this number does not depend strongly on the choice of 
this characteristic distance. Hence, the number of con- 
densed ions is reasonably well defined. Moreover, the 
excess number of multivalent counterions, which can be 
directly calculated from the experimental data, is nearly 
identical to this quantity. This point will be further ex- 
plained below. 

Fig. 8 shows the excess of 4-valent counterions 5pz{r) 
up to a distance r from the DNA axis, as a function of 
r: 

bp, (r) = 271 f r'dr' K (^') " ] (Al) 
Jo 

with the limit bpzipa) — p, of Eq. ^ The density pro- 
file was calculated using the Poisson-Boltzmann equa- 
tion, with the radius of DNA taken as c? = 10 A and 
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FIG. 8 Excess of 4-valent ions per DNA monomer, up to 
a distance r from the axis of a charged cylinder of radius 
d = lOA (modeling the DNA) as obtained using the Poisson- 
Boltzmann equation (solid line). The excess Spzif) is defined 
in Eq. lAll The number of charges per unit length on the 
cylinder is 1/a where a = 1.7 A to fit DNA values. The bulk 
densities of monovalent and multivalent ions are Ci° = 88 mM, 
cf = 0.52 mM, yielding k"^ = 10.0 A. The quantity 5p^ (solid 
line) can be compared with the total number of 4-valent ions 
(dashed line) up to a distance r from the cylinder. The dis- 
tance d + kT^ from the DNA axis is indicated by a vertical 
arrow, and characterizes the decay of the density profile far 
away from the DNA. 



with bulk densities of ions as in the last line of Table 1 : 



= c?° = 88 mM, 



0.52 mM. 



Three observations can be made. First, most but not 
all of the excess z-valent ions are localized very close to 
the DNA, at a distanc e of order A / z, wh ere A is the Gouy- 
Chapman length (see lAndelmaiJ l)l994|) '): 



A 



d 



(A2) 



where a is the average charge per unit area on the cylin- 
der surface, a = pdna/Sttc?, and pdna = 1/a is the DNA 
charge per unit length. At room temperature the Bjer- 
rum length Ib — 7 K and for DNA with 4-valent coun- 
terions A/z ~ 0.6 A. Second, the counterions within a 
layer of few times the Debye length {k^^ — 10.0 A in 
Fig. 8) neutralize the DNA charge. Nearly all the excess 
distribution is in this layer. Third, in order to estimate 
the total amount of counterions in the condensed layer 
of thickness aK~^, where a is a number of order unity, 
we need to add Spz to the bulk contribution, ira^ k~'^c^ . 
Using K from Eq. ^ the latter is equal to: 



a 

47^ 



2cf= + z{z + l)c° 



(A3) 



In experiment, is much smaller than c^j" at the onset. 



and the bulk contribution of Eq. IA3I can be neglected 
relative to pz, for a of order unity. This can be seen 



TABLE 3 Number of z-valent counterions, per 
DNA monomer, up to several different distances 
from the DNA axis, compared with ap^ 



c^mM] 


d + 10 A 


d + 20A 


d + 


d + 2k-^ 


apz 


2 


0.191 


0.193 


0.194 


0.194 


0.194 


13 


0.187 


0.190 


0.190 


0.191 


0.191 


23 


0.171 


0.172 


0.172 


0.173 


0.173 


88 


0.134 


0.135 


0.134 


0.135 


0.135 



specifically in Fig. 8 by comparing the solid and dashed 
lines. 

The outcome of the above discussion is that pz , defined 
in Eq. ^ as the excess of counterions throughout the 
cell, can be regarded, to a good approximation, as the 
total number of counterions within a condensation layer 
whose thickness is approximately the Debye length. For 
typical concentration ranges as considered here we do not 
expect that this outcome will change, even for models 
going beyond Poisson-Boltzmann theory. 

As a further demonstration, the number of multiva- 
lent counterions up to several different distances from 
the DNA is shown in Table 3, as calculated in a unit 
cell using the Poisson-Boltzmann equation. For each Cs 
in Table 1 we find the Poisson-Boltzmann density profile 
such that c'^ = Cs and pz = pi, and then calculate the 
number of multivalent ions (per DNA monomer) up to 
the following distances from the DNA radius: 10 A, 20 A, 
and 2k~^. The values of AC ; clS obtained from Eq.n 
are equal to 68, 26, 20 and 10 A for Cs = 2, 13, 23 and 
88 mM, respectively. These numbers are compared with 
ap*. All the different measures in Table 3 yield results 
that are very close to each other. 
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